Frequency-Dependent Current Noise through Quantum-Dot Spin Valves 



Matthias Braun/ Jiirgen Konig,^ and Jan Martinek^'^'^ 

^Institut fiir Theoretische Physik III, Ruhr-Universitdt Bochum, 44^80 Bochum, Germany 
^Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 
^Institut fiir Theoretische Festkorperphysik, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 
^Institute of Molecular Physics, Polish Academy of Science, 60-179 Poznan, Poland 

(Dated: February 5, 2008) 

We study frequency-dependent current noise through a single-level quantum dot connected to 
ferromagnetic leads with non-collinear magnetization. We propose to use the frequency-dependent 
Fano factor as a tool to detect single-spin dynamics in the quantum dot. Spin precession due to an 
external magnetic and/or a many-body exchange field affects the Fano factor of the system in two 
ways. First, the tendency towards spin-selective bunching of the transmitted electrons is suppressed, 
which gives rise to a reduction of the low-frequency noise. Second, the noise spectrum displays a 
resonance at the Larmor frequency, whose lineshape depends on the relative angle of the leads' 
magnetizations. 

PACS numbers: 72.25.-b, 85.75.-d, 72.70.+m, 73.63.Kv, 73.23.Hk 



I. INTRODUCTION 

The measurement of current noise reveals additional 
information about mesoscopic conductors that is not con- 
tained in the average current Current noise through 
quantum dots exposes the strongly-correlated character 
of charge transport due to Coulomb interaction, giving 
rise to phenomena such as positive cross correlations,^ 
and sub- or super-Poissonian Fano factors This is 
one motivation for the extensive theoreticalSiL2i2iiS and 
experimentaUiiiSii^ study of zero- and finite-frequency 
noise of the current through quantum dots. Furthermore, 
the finite- frequency noise provides a direct access to the 
internal dynamics of the system such as coherent oscilla- 
tions in double-dot structures ^^ii^iiSiii quantum-shuttle 
resonances^ transport through a dot with a precessing 
magnetic moment^ or back action of a detector to the 
systemi22iSii2S 

In this paper, we investigate the transport through 
a single-level quantum dot connected to ferromagnetic 
leads with non-collinear magnetizations in the limit of 
weak dot-lead coupling, see Fig. ^ Recent experimental 
approaches to contact a quantum dot to ferromagnetic 
leads involve metallic islands granular systems 
carbon nanotubes^SiSS as well as single molecules^S or 
self-assembled quantum dots^Si^ Quantum-dot spin- 
valve structures are interesting, since the presence of 
both a finite spin polarization in the leads and an ap- 
plied bias voltage induces, for a non-parallel alignment 
of the lead magnetization directions, an non-equilibrium 
spin on the quantum dot. The magnitude and direction 
of the quantum-dot spin is determined by the interplay 
of two processes: non-equilibrium spin accumulation due 
to spin injection from the leads, and spin precession due 
to an exchange field generated by the tunnel coupling 
to spin-polarized leads^^ or due to an externally applied 
magnetic field. The resulting average quantum-dot spin 
affects the dc— conductance of the device. 

While the time-averaged current is sensitive to the 
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FIG. 1: A quantum dot contacted by ferromagnetic leads with 
non-collinear magnetizations. Electrons polarized along the 
source (left) lead enter the dot. During their stay on the dot, 
the spins precess in the many-body exchange field -Bl + -Br, 
which arises from the tunnel coupling to the left and right 
lead, and an applied magnetic field -Bapp. Due to magne- 
toresistance effects this precession modulates the tunnel-out 
probability to the drain (right) lead, giving rise to a signal in 
the power spectrum of the current noise. 



time-averaged dot spin, the time-resolved dynamics of 
the dot spin is provided by the power spectrum of the 
current noise. It will show a signature at the frequency 
that is associated with the precession of the quantum- 
dot spin due to the sum of exchange and external mag- 
netic field. This can be understood by looking at the 
tunneling-out current to the drain (right) lead as a func- 
tion of the time after the quantum-dot electron had tun- 
neled in from the source (left) lead. The spin of the 
incoming electron, defined by the source-lead magnetiza- 
tion direction, precesses about the sum of exchange and 
external magnetic field as long as it stays in the dot. 
Since the tunneling-out rate depends on the relative ori- 
entation of the quantum-dot spin to the drain-lead mag- 
netization direction, the spin precession leads to a pe- 
riodic oscillation of the tunneling-out probability. The 
period of the oscillation is defined by the inverse pre- 
cession frequency, and the phase is given by the relative 
orientation of the source- and drain-lead magnetization 
direction. As a consequence, the signature in the power 
spectrum of the current noise at the Larmor frequency 
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gradually changes from a peak to a dip as a function of 
angle between source- and drain-lead magnetization. 

Also the zero- frequency part of the current-noise power 
spectrum is affected by the internal dynamics of the 
quantum-dot spin. By coupling a quantum dot to spin- 
polarized electrodes, the dwell time of the electrons in 
the dot becomes spin dependent. It is knowni^ that 
this spin dependence of the dwell times yields a bunch- 
ing of the transferred electrons, which leads to an increase 
of the shot noise. A precession of the quantum-dot spin 
due to exchange and external magnetic field weakens the 
tendency towards bunching, leading to a reduction of the 
low-frequency noise. 

The aim of this paper is to perform a systematic study 
of the frequency-dependent current noise of a quantum- 
dot spin valve in the limit of weak dot-lead coupling 
in order to illustrate the effects formulated above. In 
Sec.^we define the model of a quantum-dot spin valve, 
as shown in Fig. ^ Sec. IIIIl we extend a previously- 
developed diagrammatic real-time technique^i to eval- 
uate frequency-dependent current noise, as it has been 
similarly done for metallic (non-magnetic) single-electron 
transistors ^2i2& The results for the quantum-dot spin 
valve are discussed in Sec. IIVI followed by the Conclu- 
sions in Sec. 

II. MODEL SYSTEM 

The Hamiltonian for the quantum-dot spin valve, i.e., 
a quantum dot coupled to ferromagnetic leads, is given 
by the sum 

H = Hi^^Hn^Hj,^HT. (2.1) 

The single-level quantum dot is modeled by an Anderson 
impurity, 

Hd= + U n^n^ , (2.2) 

where cj. and are the fermion creation and annihila- 
tion operators of the dot electrons, and Ua- = c'lc^. The 
single-particle level at the energy e, measured relative to 
the equilibrium Fermi energy of the leads, may be split 
due to an external magnetic field, = e -\- A/2 and 
£| = e — A/2 with Zeeman energy A = ^/ie^ext- Double 
occupancy of the dot costs the charging energy U /cbT. 

The ferromagnetic leads (r = L/R) are treated as 
reservoirs of non-interacting fermions, 

/c,q;=± 

By choosing the quantization axis of each lead parallel 
to their direction of magnetization fi^, the property of 
ferromagnetism can be included by assuming different 
density of states ^q, for majority {a = +) and minority 
{a = —) electrons. An applied bias voltage is incorpo- 
rated by a symmetric shift of the chemical potential by 



Ml/r = ±6^/2 in the left and right lead, which enter the 
Fermi functions fr{E) = f{E — jUr)- 

The magnetization directions of the left and right lead 
and the external magnetic field are, in general, non- 
collinear, i.e., in the Hamiltonians for the three subsys- 
tems we have chosen different spin quantization axes. 
To describe spin-conserving tunneling, one must include 
SU{2) rotation matrices U^^ in the tunneling Hamilto- 
nian 

Ht= traU^K^c^ + H.c. . (2.4) 

r,k,aa 

For simplicity we use leads with energy-independent den- 
sity of states and barriers with energy-independent 
tunnel amplitudes tr. With these assumptions, the de- 
gree of lead polarization p = (^+ —£,-)/ (^+ + ^- ) as well 
as the coupling constants F^ = Xla^i 27r|tr/\/2p^^a do 
not depend on energy. 

III. DIAGRAMMATIC TECHNIQUE 

The dynamics of the quantum-dot spin valve is deter- 
mined by the time evolution of the total density matrix. 
Since the leads are modeled by non-interacting fermions, 
which always stay in equilibrium, we can integrate out 
the degrees of freedom in the leads, and only need to 
consider the time evolution of the reduced density matrix 
p{t) of the quantum dot, which contains the information 
about both the charge and spin state of latter. In the 
following three subsections, we formulate the derivation 
for the stationary density matrix, the dc— current and the 
finite-frequency current-current correlation function. Af- 
terwards, in Sec. IIII Dl we specify the obtained formulas 
for the limit of weak dot-lead coupling, i.e. we perfrom 
a systematic lowest-order perturbation expansion in the 
tunnel coupling strength F = Fl + Fr. 

A. Density matrix 

The quantum-statistical average of the charge and spin 
on the quantum dot at time t is encoded in the reduced 
density matrix p{t). Its time evolution is governed by 
the propagator n(t,to), 

p{t)=U{t,to) •p(to). (3.1) 

Since p is a matrix, the propagator 11) must be a ten- 
sor of rank four. A diagrammatic representation of this 
equation (see also Ref. 35) is depicted in Fig. [3 The 
upper/lower horizontal line represents the propagation 
of the individual dot states forward/backward in (real) 
time, i.e. along a Keldysh time contour tx- 

In order to find the stationary density matrix for 
a system, which is described by a time-independent 
Hamiltonian, we consider the limit to — oo. There 
is some characteristic time after which the system 
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FIG. 2: The density matrix evolves in time with the propa- 
gator n, which is a tensor of rank four. 



loses the information about its initial density matrix 
Pini = limto^-oo p(^o)- We can, therefore, choose with- 
out loss of generality (pini)^^^ = ^xi,xo^r7i,xo with an 
arbitrarily-picked state xo, to get for the stationary (non- 
equilibrium) density matrix 
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independent of xo- Here, for time-translation invariant 
systems, the propagator 11 (t, to) depends only on the dif- 
ference of the time arguments {t — to). For the following, 
it is convenient to express the propagator in frequency 
representation n(a;) = J^^dtll{—t) ex.p[i{Lj — 
iO~^)t]. It can be constructed by the Dyson equation 

= [Uo-\u;)-W{u;)]-' . (3.3) 

The full propagator 11 (cj) depends on the free propaga- 
tor Ilo{u;) and the irreducible self energies W{u;)^ which 
describes the influence of tunneling events between the 
dot and the leads. The Dyson equation is diagrammat- 
ically represented in Fig. |31 The frequency argument of 
the Laplace transformation appears in this diagrammatic 
languageSi as additional horizontal bosonic line trans- 
porting energy hco. 
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FIG. 3: Diagrammatical representation of the Dyson equa- 
tion for the propagator. The self energy W sums up all ir- 
reducible tunnel diagrams. With VF(cj), we label the self 
energy, together with the parallel running frequency line cj. 



leads are in equilibrium, their non-interacting fermionic 
degrees of freedom can be integrated out. Thereby two 
tunnel Hamiltonians each get contracted, symbolized by 
a line. Each line is associated with one tunnel event, 
transferring one particle and a frequency /energy from 
one vertex to the other. Therefore the lines have a de- 
fined direction and bear one order of the coupling con- 
stant r = Fl + Fr. We define the self energy W{uj) 
as the sum of all irreducible tunnel diagrams (diagrams, 
which can not be cut at any real time, i.e. cut vertically, 
without cutting one tunneling line). 

In Sec, nil D1 we will then restrict our otherwise general 
calculation to the lowest-order expansion in F, i.e. we will 
include only diagrams with one tunnel line \ViW(uj). A 
detailed description of how to calculate these lowest-order 
self energies as well as example calculations of W for the 
system under consideration can be found in Ref. 

To solve for the stationary density matrix pst, we 
rewrite the Dyson Eq. as (no(cj)"^ - Vr(cj))n(cj) = 
1, multiply both sides of the equation with uj^ use 
the final value theorem lim^^o(^^^ + 0+)n(c(;) = 
lim^^oo n(t), similar as for Laplace transformations, and 
employ Eq. (|3.2^) . to get the generalized master equation 
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(3.5) 



together with the normalization condition Tr[pst] = 1. 

The structure of Eq. (|3.5|) motivates the interpreta- 
tion of the self energy W{ix) = 0) as generalized transi- 
tion rates. However, the self energy does not only de- 
scribe real particle transfer between leads and dot, but 
it also accounts for tunneling-induced renormalization 
effects. It was shown in Refs. ^,37,38.39, that these 
level renormalization effects may affect even the lowest- 
order contribution to the conductance. Therefore, a ne- 
glect of these renormalizations would break the consis- 
tancy of the lowest-order expansion in the tunnel cou- 
pling strength. MiiSiii Recently, the frequency-dependent 
current noise of a quantum-dot spin valve structure was 
discussed in Ref. [42, in the limit of infinite bias voltage, 
where these level renormalizations can be neglected. One 
of the main advantages of the approach presented here 
is, that a rigid systematic computation of the generalized 
transition rates is possible, which include all renormaliza- 
tion effects. Therefore our approach is valid for arbitrary 
bias voltages. 



The free propagator (without tunneling) is given by 



Current 



(3.4) 



where (e^) is the energy of the dot state x iv)- Tun- 
neling between the dot and the leads introduce the irre- 
ducible self energies W{uj). We calculate W{uj) in a per- 
turbation expansion in the tunnel Hamiltonian Eq. (|2.4|) . 
Each tunnel Hamiltonian generates one vertex (filled cir- 
cle), on the Keldysh time contour tx, see Fig. 01 Since the 



The current through barrier r = L, R is defined as the 

change of charge erir = e ^Ika^rka ^^^^ ^ 
tunneling, described by the operator 



drir e 



(3.6) 



We define the operator for the current through the dot 
as / = (/l — /r)/2. Each term of the resulting current 
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operator does contain a product of a lead and a dot op- 
erator. By integrating out the lead degrees of freedom, 
the current vertex (open circle) gets connected to a tun- 
nel vertex by a contraction line as depicted in Fig. ^ 
Thereby the tunnel vertex can be either on the upper or 
lower time contour line. 

!={!)= Tr[/pst] = (pstT' W i 
\ y m T]-^ ^ 

FIG. 4: Diagrammatic representation of the current. By 
integrating out the lead degrees of freedom, the current vertex 
(open circle) gets contracted to one of the tunnel vertex in a 
self energy W(io — 0). 

To present a systematic way to calculate the current, 
we can utilize the close similarity of the tunnel Hamilto- 
nian in Eq. (|2.4|) and the current operator in Eq. (|3.6|) . 
Both differ only by the prefactor e/h and possibly by 
additional minus signs. 

Following the work of Thielmann et al. ^ we define the 
object W^^^^^ as the sums of all possible realizations of 
replacing one tunnel vertex (filled circle) by a current 
vertex (open circle) in the self energy W^_^^^ compare 
Fig. [51 In technical terms, this means that each diagram 
is multiplied by a prefactor, determined by the position 
of the current vertex inside the diagram. If the current 
vertex is on the upper (lower) Keldysh time branch, and 
describes a particle tunneling into the right (left) lead 
or out of the left lead (right), multiply the diagram by 
+ 1/2, otherwise by —1/2. For clarity, we keep the factor 
e/h separate. For the detailed technical procedure of the 
replacement as well as the rules to construct and calculate 
the self energies, we refer to Ref. 0. The average of one 
current operator, i.e. the dc— current fiowing through the 
system is then given by 

I={i) = -^l^[W\LO = Q)p,,]. (3.7) 

The trace selects the diagonal matrix elements, which 
regards that the Keldysh line must be closed at the end 
of the diagram, see Fig.[5J requiring that the dot state of 
the upper and lower time branch match. 




FIG. 5: Reformulation of the current as function oiW^ {uj = 
0), the self energy with one tunnel vertex replaced by a current 
vertex. 

To see, that the diagrams in Fig. ^ and Fig. O are 
equal, one must consider, that all diagrams, where the 
rightmost vertex is a tunnel vertex will cancel each other 



when performing the trace. This happens, since by mov- 
ing the rightmost tunnel vertex from the upper (lower) to 
the lower (upper) Keldysh time line, the diagram aquires 
only a minus sign^ 



C. Current-current correlation 

We define the frequency-dependent noise as the Fourier 
transform oi S{t) = (/(t)/(0)) + (/(0)/(t)) - 2(/)2, which 
can be written as 

S{uo) = / dt((/(t)/(0)) + (/(0)/(t)))(e-^"^ + e+^"^) 

. (3.8) 

We restrict our discussion to the above defined sym- 
metrized and, therefore, real noise since it can be mea- 
sured by a classical detector. The unsymmetrized noise 
would have an additional complex component, describing 
absorption and emission processes,^ that depend on the 
specifics of the detector. 

Since the dot-lead interface capacitances are much less 
sensitive to the contact geometry than the tunnel cou- 
plings Fl/j^, we assume an equal capacitance of the 
left and right interface, while still allowing for differ- 
ent tunnel-coupling strengths. Following the Ramon- 
Shockley^^ theorem, we have then to define current oper- 
ator also symmetrized with respect to the left and right 
interface as already done in Sec. IIII Bl 

The diagrammatic calculation of the current-current 
correlation function is now straightforward. Instead of 
replacing one tunnel vertex by a current vertex on the 
Keldysh time contour, as for the average current, one 
must replace two vertices. The additional frequency uo 
of the Fourier transformation in Eq. (|3.8|) can be incor- 
porated in the diagrams as an additional bosonic energy 
line (dashed) running from t to 0, i.e. between the two 
current vertices. This line must not be confused with a 
tunnel line, since it only transfers energy hjj^ and no par- 
ticle. By introducing the self energy W all diagrams of 
the current-current correlation function can be grouped 
in two different classes^^^^ as shown in Fig. [SI Either 
both current vertices are incorporated in the same irre- 
ducible block diagram, or into two different ones that are 
separated by the propagator Yi{uo). 

[mm) + (i{o)i{t)))e-'-' = 

^'■^-- ui 

= p,t w j +Pst ~^w n 1^"^ + 

^ f^- d 

+ diagrams with one or both vertices on the lower line 

FIG. 6: Regrouping of the noise expansion by introducing 
the irreducible self energy and the propagator 11(0;). 
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The order of the current operator on the Keldysh con- 
tour is determined by its ordering in the correlator, so the 
current operator at time hes on the upper branch for 
(/(t)/(0)) and on the lower branch for (/(0)/(t)). Since in 
Eq. dSSI) 

we defined the noise symmetrized with respect 
to the operator ordering, we just allow every combination 
of current vertex replacements in the W^s. This includes 
also diagrams where one or both vertices are located on 
the lower time contour (this type of diagrams are not 
explicitly drawn in Fig.|H|). 

By including the current vertices and the frequency 
line in the self energies, three variants of the self energy 
W are generated. The objects W^{u) and W^{ij) are 
the sum of all irreducible diagrams, where one tunnel 
vertex is replaced by a current vertex in any topological 
different way. The subindex > (<) indicates, that the 
frequency line connected to the current vertex leaves or 
enters the diagram to the right (left) side. In the zero- 
frequency limit, the two objects become equal Wy{uj = 
0) = W^ico = 0) = W^. 

The third object W^^{u;) sums irreducible diagrams 
with two tunnel vertices replaced each by a current ver- 
tex in any topological different way. The current vertices 
are connected by the frequency line co. The diagrammat- 
ical picture of the objects W{uj), W^^{uj), W>(cc;), and 
W^{uj) are shown in Fig. [3 
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FIG. 7: Different variations of the self energy W . 

With these definitions the diagrams for the frequency- 
dependent noise in Fig.lHlcan be directly translated into 
the formula 



S{l0) = —TT[W"iL0)pst + Wi{L0)UiL0)Wl{L0)pst] 



-2Tr5iLu) {If + (w ^ -uj) . 



(3.9) 



We remark that the first line in Eq. (|3.9|) diverges as Cc; ^ 
0. While the VF's are regular for ^ 0, the propagator 
n(cj) goes as i/{—hLJ-\- iO~^) times Il{t ^ oo), which 
is related to Pst via Eq. (|3.2^) . In the limit uj ^ the 



propagator therefore yields both a delta function 5{oo) 
and a l/o; divergence. For the full expression of the noise, 
these divergences are canceled by the delta-function term 
in the second line of Eq. and by the terms with 

UJ —uj^ respectively. As a consequence, S{uj) remains 
regular also in the limit ^ 0. 



D. Low- frequency noise in the sequential-tunnel 
limit 

Equation (|3.9|) is the general expression for the 
frequency-dependent current noise. In the following pa- 
per, we consider only the limit of weak dot-lead tunnel 
coupling, r <C /cbT, and therefore include only diagrams 
with at most one tunnel line in the VF's. However, this 
procedure is not a consistent expansion scheme for the 
noise S{uj) itself. By expanding the W's up to linear 
order in F, the result of Eq. (|3.9|) is the consistent noise 
linear in F plus some higher-order contributions propor- 
tional to F^. Since co-tunnel processes also give rise to 
quadratic contributions, we have to discard these terms 
as long as we neglect the quadratic cotunnel contribu- 
tions ofW. If one is interested in the noise up to second 
order in F, then these higher-order terms generated by 
lower-order VF's are of course an essential part of the 
result."^^ 

Further, we are looking for signatures of the inter- 
nal charge and spin dynamics of the quantum-dot in the 
frequency-dependent current noise. Therefore - if we ne- 
glect external magnetic fields at this point - we concen- 
trate on frequencies that are at most of the same order 
of the tunnel coupling F. If we limit the range in which 
we want to calculate the current noise to ficj < F, we 
can neglect the frequency dependence of the W's. Each 
correction of the W's would scale at least with ujT ^ F^, 
making them as important as the neglected co-tunnel 
processes. 

The neglect of the terms in W which are at least linear 
in frequency has two main advantages. First, it consid- 
erably simplifies the calculation of the W's. Second, it 
automatically removes the quadratic parts of the noise, 
so Eq. (|3.9|) gives a result consistent in linear order in F. 
In this low- frequency limit, the noise can then be written 



-27r5{uj) {If + (cj ^ -uj) 



(3.10) 



where = Wl{uj = 0) = W^{uj = 0), W = W{uj = 
0), and W^^ W^^{uj 0). This means, that the 
bosonic frequency lines uj in the diagrams as shown in 
Fig. El can be neglected. The only remaining frequency- 
dependent part is the free propagator Tlo{uj). 

This formalism, of course, reproduces the noise spec- 
trum of a single-level quantum dot connected to normal 
leads as known from literature i If one can approximate 
the Fermi functions by one or zero only, i.e. if the dot 
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levels are away from the Fermi edges of the leads the 
Fano factor shows a Lorentzian dependence on the noise 
frequency uo 



F{lo) 



(2rL - rR)2 



S{lo) _ 1 

2ei ~ 2 + (2rL + rR)2 + {hwY 



(3.11) 



for a bias voltage allowing only an empty or singly- 
occupied dot, and 



F{lo) 



1 



2 (Fl + rR)2 + (;k^)2 



(3.12) 



for higher bias voltages, when double occupation is also 
allowed. 



E. Technical summary 

The technical scheme for calculating the zero- and low- 
frequency current noise is the following: First the objects 
W ^ and W^^ must be calculated in the a; = limit, 
using the diagrammatic approach, see Ref. 32. 

In the next step, we calculate the reduced density ma- 
trix p of a single- level quantum dot, which is a 4 x 4 
matrix, 



/ \ 

pj 

p| p| 

V pd / 



(3.13) 



since the dot can be either empty (x = O)^ occupied with 
a spin- up (x =T) oi" spin-down (x =i) electron, or 
doubly occupied (x = d). The diagonal elements of the 
matrix can be interpreted as the probability to find the 
dot in the respective state, while the inner 2x2 matrix 
is the SU{2) representation of the average spin on the 
dot. All off-diagonal elements connecting different charge 
states are prohibited by charge conservation. 

For technical reasons it is convenient, to ex- 
press the density matrix as vector: pst = 
(Po^^/,^/,^d,^/,^/)^. Then the forth-order 
tensors W^s and n(cj)'s are only 6x6 matrices, see 
App. ^ and standard computer implemented matrix 
operations can be used. It is worth to point out, that in 
the vector notation, the trace for example in Eq. (|3.10|) 
is then not the sum of all elements of the resulting 
vector as assumed by Ref. [lEl but only the sum of 
the first four entries. These elements correspond to 
the diagonal entries of the final density matrix. In the 
notation of Ref. 0, this can be achieved by the vector 
e^ = (1,1,1,1,0,0). 

The stationary density matrix follows from the master 
Eq. est = -i(e,-ex)(P^t)^+E;,„,, W>^^^^ {pst)f, un- 
der the constraint of probability normalization • pst = 
1. The average dc— current through the system is given 
by / = e/{2h)e^ • • pst- In the low-frequency limit 



the frequency-dependent propagator n(cj) can be con- 
structed from the frequency-dependent free propagator 
no(a;) and the fieqnency- independent self energy W{u; = 
0). The low-frequency noise is then given by the matrix 
multiplication = e'^ / {2h)e^ -{W^^ + W^U{ij)W^)- 
Pst ^ — where the zO+ in the denominator of the 
propagator is dropped, since the term arising from the 
zO+ contribution cancels the delta function in Eq. (|3.9^) . 



IV. RESULTS 

In this section, we discuss our results for zero- and 
finite-frequency current noise in a quantum dot con- 
nected to ferromagnetic leads with non-collinear magne- 
tizations. The relative energies of a single-level dot is 
sketched in Fig. [HI 
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FIG. 8: Sketch of different energies involved. Since we assume 
equal tunnel interface capacities, the voltage drop on the left 
and right side is symmetric. 



We always assume /cbT ^ F, and that the single- 
particle state is above the equilibrium Fermi energy of 
the leads, otherwise higher-order tunnel processes could 
become important ^i^^i^Siii 



A. Zero-frequency noise 

We start our discussion with the zero-frequency noise. 
In Fig. El we plot results for F{uj ^ 0) = S{uj 
0)/(2e/), i.e. the zero-frequency Fano factor for the 
quantum dot contacted by ferromagnetic leads. In 
Fig. Et), the leads are aligned parallel. For eV/2 < £, 
when the dot level is above the lead Fermi energies, the 
dot is predominantly empty, and interaction effects are 
negligible leading to a Fano factor of 1. In the volt- 
age window e < eV/2 < e -\- when the dot can 
only be empty or singly occupied, we can observe super- 
Poissonian noise due to dynamical spin blockad^iSi24 for 
sufficiently high lead polarization. The minority spins 
have a much longer dwell time inside the dot than the 
majority spins. In this way, they effectively chop the 
current leading to bunches of majority spins. While the 
current in this regime / = 2rLrR/(2rL + Fr) does not 
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depend on the polarization p of the leads, the Fano factor 

4 i±4r2 + 

F{0) = " (4.1) 

(2rL + Tr)^ 

even diverges for p ^ 1. If the voltage exceeds the 
value necessary to occupy the dot with two electrons 
{eV/2 > £ + the noise is no longer sensitive to a 
lead polarization. 




eV / k^T 



FIG. 9: Zero- frequency current noise through a quantum dot 
spin valve. In panel a), the lead magnetizations are aligned 
parallel, in panel b) anti-parallel, and in panel c) the lead mag- 
netizations enclose an angle of 7r/2. The different lines cor- 
respond to different values of the lead polarization p. Other 
parameters are e = 10/cbT, U = SO/cbT, and Fl = 2Fr 



lead and weakly coupling to the drain lead, therefore the 
Fano factor approaches unity. 

If the leads are non-collinearly aligned, for example 
enclose an angle (j) = 7r/2 in Fig. [Uh), a qualitatively 
different behavior can be observed. Now, the typical 
Coulomb plateaus are modulated. This shape arises, 
since the dot spin starts to precess around the lead mag- 
netizations. The tunnel coupling between the ferromag- 
netic lead r =L/R and the dot induces the exchange field 
contribution22i2^ 



^ J 



duj 



1 - frico) 



u 



uo ■ 



, (4.2) 



generating an intrinsic spin precession of the dot spin 
around the lead magnetizations. This exchange field au- 
tomatically appears in a rigid calculation of the general- 
ized transition rates W. 

The intrinsic spin precession due to the exchange field 
counteracts the dynamical spin blockade. The exchange 
coupling to one lead is maximal, if its Fermi energy co- 
incides with the dot energy levels, i.e. the coupling to 
the source lead is maximal at the voltages eV/2 = e and 
eV/2 = e-\-U and changes its sign in between. Therefore 
the reduction of the Fano factor is non-monotonic, and 
so is the variation of the Coulomb plateaus. It is worth 
to point out, that to observe this spin precession mech- 
anism in the conductance of the device a relatively high 
spin polarization of the leads is required. But the noise is 
much more sensitive to this effect than the conductance, 
that a polarization as expected for Fe, Co, or Ni4& is well 
sufficient. 

The zero-frequency Fano factor as a function of the an- 
gle (j) between the two lead magnetization vectors is plot- 
ted in Fig. The black lines are for the bias voltage 
eV = 50/cbT, where the exchange field influence is weak, 
while the gray lines is for the bias voltage eV = SOkBT. 
Since both voltages are within the voltage window allow- 
ing only single occupation of the dot, compare Fig. El the 
tunnel rates do not change significantly within this volt- 
age range. Only the exchange field varies with voltage. 
Since the exchange field suppresses bunching due to spin 
precession, the black and gray curves split. 

For = and (j) = tt the accumulated spin is 
collinearly aligned with the exchange field, and no spin 
precession arises. 



Also in the case of anti-parallel aligned leads, the Fano 
factor rises in the voltage regime e < eV/2 < e -\- U as 
seen in Fig. [Ub). The dot is primarily occupied with 
an electron with majority spin of the source lead, i.e. 
minority spin for the drain lead, since this spin has the 
longest dwell time. If the electron tunnels to the drain 
lead, it gets predominantly replaced by a majority spin 
of the source lead. For a high enough lead polarization, 
only one spin component becomes important. Further 
this spin component is strongly coupled to the source 



B. Finite- frequency noise and weak magnetic fields 

The dc— conductance of the quantum-dot spin valve is 
a direct measure of the time- aver aged spin in the dot. On 
the other side, the power spectrum of the current noise 
can also measure the time-dependent dynamics of the 
individual electron spins in the dot. The spin precesses in 
the exchange field as well as an external magnetic field. 
This gives rise to a signal in the frequency-dependent 
noise at the Larmor frequency of the total field. 
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FIG. 10: Fano factor of a quantum dot spin valve as a 
function of the angle 0, enclosed by the lead magnetiza- 
tions. The black lines are for a bias voltage eV = 50/cbT, 
where the exchange field is very weak, while the gray curves 
are for eV = 30/cbT, where the exchange field is more pro- 
nounced. Further parameters are s — 10/cbT, U = SO/crT, 
and Fl = 2Fr 



(zeroth-order) stationary density matrix pst, via the mas- 
ter Eq. (ESI). 

The numerical results are plotted in Fig. and 
Fig. El In Fig. the magnetizations of the leads are 
aligned parallel, and a magnetic field is applied perpen- 
dicular to the lead magnetizations. With parallel aligned 
leads and equal polarizations in both leads, no average 
spin accumulates on the dot, and therefore the current- 
voltage characteristic as shown in the inset of Fig. ITTl 
shows neither magnetoresistance nor the Hanle effect if a 
transverse magnetic field is applied. In contrast to the 
conductance, which depends on the average dot spin only, 
the frequency-dependent noise is sensitive to the time- 
dependent dynamics of the spin on the dot. Therefore 
the field-induced spin precession is visible in the noise 
power spectrum. For B = the Fano factor shows a 
Lorentzian dependence of the noise frequency. Thereby 
the Fano factor exceeds unity due to the bunching effect, 
as discussed in Sec. IIV Al 

With increasing magnetic field, spin precession lifts the 
dynamical spin blockade inside the dot, and the Fano 
factor decreases at ~ 0. 



By including an external magnetic field in the noise 
calculation, one has to distinguish two different parame- 
ter regimes: either the Zeeman splitting A = ^/iB^ext is 
of the same order of magnitude as the level broadening 
A rL,rR, or it significantly exceeds the tunnel cou- 
pling A rL,rR. In this section we focus on the first 
case, while the latter case is treated in Sec. IIV (^1 

By choosing the spin-quantization axis of the dot sub- 
system parallel to the external magnetic field, the mag- 
netic field only induces a Zeeman splitting of the single- 
particle level e in = e -\- A/2 and ei = e — A/2. Since 
A Fl, Fr, we can expand the Ws also in A and keep 
only the zeroth-order terms, since each correction of the 
self energies would be proportional to A • F F^. The 
Zeeman splitting must only be considered for the free 
propagator. With Eq. (|3.4|) . the propagator is then given 
by 



no(c.^) 



f CO 

uj 



CO 

cc; 

V 



(4.3) 



A/ 



where we already dropped the +zO+ in the denominator, 
and use the matrix notation as introduced in Sec. IIII El 
The two last rows of this matrix govern the time evo- 



lution of p| and p|, representing the spin components 



transverse to the quantization axis, i.e. transverse to the 
applied magnetic field. The change of the denominator 
by the Zeeman energy A describes just the precession 
movement of the transverse spin component. Since the 
free propagator no(a;) is a function of A, the Zeeman 
energy modifies the full propagator Il{co) as well as the 




co/F 



FIG. 11: Frequency-dependent Fano factor of a quantum dot 
connected to parallel aligned leads for various perpendicular 
applied external magnetic fields. The parameters are p = 0.5, 
e = 10/cbT, U = 30A;bT, eV = 40A:bT and Fl = 2rR. The 
inset shows the current bias- volt age characteristic, which does 
not depend on the applied magnetic field. 



Further a resonance line evolves approximatively at the 
Larmor frequency of the applied magnetic field. The line 
width of the resonance is given by the damping due to 
tunnel events. If the dot can only be singly occupied, the 
damping coefficient equals the tunnel-out rate Fr, if the 
dot can also be doubly occupied, also tunnel-in events 
contribute. 

The deviation of the resonance line position from the 
Larmor frequency, one would expect by considering the 
applied magnetic field only, is caused by the exchange in- 
teraction. The spin inside the dot precesses in the total 
field containing the external magnetic field and the ex- 
change field. Dependent on their relative orientation. 
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the exchange field can increase or decrease the total field 
strength. 

Since the exchange field is a function of the applied 
bias voltage, the resonance peak is shifted by changing 
the bias voltage. In Fig. the finite frequency noise 
for a quantum dot is plotted, where the lead's magne- 
tizations enclose an angle (j) = 7r/2, i.e. their magne- 
tizations are perpendicular to each other. Further an 
external magnetic field is applied parallel to the source 
lead magnetization. The exchange fields originating from 
the leads has to be added to the external field. By vary- 
ing the bias voltage (without significantly changing the 
transition rates, as indicated in the inset of Fig. [T2j) the 
exchange field varies, and the position of the resonance 
peak is shifted. 



0.8 



0.75- 



^ 0.7 
in 



0.65- 



0.6 




0.5 



co/F 



1.5 



FIG. 12: The Fano factor of a quantum dot spin valve as 
a function of the noise frequency. The lead magnetizations 
enclose an angle = 7r/2 and an external magnetic field 
^/XB^ext = l/2r is applied parallel to the source lead magne- 
tization. The vertical gray line marks the Larmor frequency 
given by the external magnetic field only. For the three differ- 
ent bias voltages, eV = 30/cbT (dot-dot-dashed), eV = 45A:bT 
(dashed) and eV = 70/cbT (solid), the strength of the ex- 
change field varies, and so does the position of the resonance 
peak. Other system parameters are p = 0.5, e = 10/cbT, 
U = SO/cbT, and Fl = 2Fr 



C. Limit of strong magnetic fields 

In this section, we discuss the case of an applied mag- 
netic field, where the Zeeman energy A = ^/iB^ext ^ 
Tl, Tr exceeds the tunnel coupling strength. As a simpli- 
fication, we can consider the tunnel rates (i.e. the VF's) 
still as independent of A as well as of c^;. This assump- 
tion is justified, if the distance between the quantum dot 
states and lead Fermi surfaces well exceeds temperature 
/cbT, the Zeeman splitting A and the noise frequency huo. 

For a clear analytic expressions, we expand the sta- 
tionary density matrix in zeroth order in F/A. Further 
we consider only the noise frequency range co = A ± F. 
In this regime the first five diagonal entries of the free 



propagator in Eq. (|4.3|) can be treated as zeroth order in 
F, i.e. their contribution drops out for the lowest-order 
noise, only the last entry 1/ {uo — A) ^ l/T is kept. This 
considerably simplifies the calculation, since all bunching 
effects and the exchange field components perpendicular 
to the external field can be neglected. 

Let us consider a single-level quantum dot with such 
an applied voltage, where approximately /l(^) = 1 and 
/l(£: + /7) = /r(£) = fR{e + [/) = 0, i.e. the applied bias 
voltage allows only an empty or singly-occupied dot. For 
an external applied magnetic field perpendicular to both 
lead magnetizations ^ext 't^L^'i^R the Fano factor 



cos ( 



rR(^ — A) sin( 



(4.4) 



shows a resonance signal at the Larmor frequency hco = 
A. By representing the frequency-dependent Fano factor 
as an integral over time. 



F{lo) 



Ah 



5R 



L 



oo 

dt e~ 



rn/h-t -i{A/h-t- 



,(4.5) 



the discussion of the functional form becomes more trans- 
parent. At t = an electron tunnels from the source 
(left) lead in the dot. This electron decays on average 
to the drain (right) lead on the time scale ^i/Fr. During 
its dwell time the electron precesses inside the dot with 
the Larmor frequency A/h. This precession modulates 
the decay rate, due to magnetoresistance effects. The 
tunnel-out event is more likely, if the spin is aligned par- 
allel to the drain lead magnetization than if anti-parallel 
aligned. The phase of this modulation is given by the 
relative angle of the lead magnetizations, and the effect 
can give rise to an absorption or dispersion line shape, 
see Fig. [HI 



(N 



0.8 
0.7 
0.6 

;o.5 

0.4 
0.3 



_ 1 1 1 1 


1 


1 H A B®Tlp)^t: 


/ \ FM 


QD FM _ 






^^^^ \ / 






(t> = 




- - (|) = 71/2 


1 T 1 1 


(^ = n 



5 
(co-A)/Fr 



10 



FIG. 13: Fano factor as a function of noise frequency for 
different angles (j) in the frequency range of the Larmor fre- 
quency. The applied voltage does only allow single occupation 
of the dot. Other system parameters are as in Fig. 1121 

By shifting the gate voltage such that /l(^) = /l(^ + 
U) = fR{e) = 1 and /r(£H-[/) = 0, the dot will always be 
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at least occupied by one electron. Then the noise shows 
the same resonance, only Fr and (j) must be replaced by 
Fl and —(j). 

If the leads are aligned parallel, the electron will leave 
the dot primary directly after the tunnel-in event, or af- 
ter one revolution, i.e. the decay is modulated with a 
cosine function. If the leads are aligned perpendicular 
to each other, then the electron must be rotated by the 
angle 7r/2 (or 37r/2) before the maximum probability for 
the tunneling-out event is reached. The decay is then 
modulated by a (minus) sine function. 

The phase dependence of the noise resonance is also 
predicted for a double-dot systemJ^ii^ Let us consider 
two dots connected in series, see Fig. ^^), and an elec- 
tron from the left (source) electrode enters the left dot. 
Since this is not an eigenstate of the isolated double-dot 
system, the electron coherently oscillates between the two 
dots with the frequency cjr. After the time t = tt/cjr, 
the electron is in the right dot and can tunnel to the 
drain lead. This corresponds to the (j) = tt case resulting 
a dip in the noise. The realization of the = case 
would be a double dot where the left (source) and right 
(drain) lead is contacted to the same dot, see Fig. I14B). 
Here the electron must stay a multiple of 27r /ujr inside 
the double dot to tunnel to the drain lead, giving a peak 
in the frequency noise spectrum. Other values of (j) have 
no double-dot-system analogon. 

A. ) = TT 

5^o(r L)c5fe(R 3oc5^ 

^ = t = n/ojR 

B. )4> = 

o ® □ o 
2L^g i)B(i ggg i)0^ 

^ = t^n/uR t = 27r/a;R 

FIG. 14: The double dot analog for the decay phase shift 
= and = TT of the electrons. 



D. Influence of spin relaxation 

The density matrix approach offers a way to phe- 
nomenologically include spin relaxation by supplement- 
ing the matrix W by 



W = w 
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■ (4.6) 



The entries in the lower right corner of Eq. (|4.6|) de- 
scribe the exponential decay of the transverse spin com- 
ponents on the time scale T2, and the block in the up- 
per left corner describes an equilibration of the occu- 
pation probability for spin up and down. If one de- 
fine the average spin vector on the quantum dot by 
5' = (p| + Pp^pj — w\ip\ — the master Eq. (|3.5|) 

becomes a Bloch equation^ The new term in Eq. (|4.6|) 
introduces an additional exponential decay term in this 
Bloch equation. In the limit of weak Zeeman splitting as 
discussed throughout the paper, Ti and T2 become equal, 
and W includes an isotropic exponential damping of the 
spin on the dot. Thereby the master equation describ- 
ing the change of the probability dt{Pj + P| ) for single 
occupation is not affected by this relaxation term. 

The modified rate matrix W enters the noise calcula- 
tion via the calculation of the stationary density matrix 
and via the propagator 11 (cj). The numerical solution for 
the case of parallel aligned lead magnetizations is plotted 
in Fig. [T^ 
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FIG. 15: Frequency dependence of the Fano factor if the leads 
are aligned parallel. With increasing spin relaxation, the spin 
blockade and therefore the bunching effect is reduced. Other 
system parameters are as in Fig. ITTl 

With increasing the spin decoherence, the spin related 
effects decrease, which is the expected behavior for spin- 
decoherence. To completely suppress the spin related 
effects the spin life time must significantly exceed the 
inverse tunnel coupling, i.e. the spin related effects are 
not very fragile against spin-decoherence. 

Several articlesiSiiSiiS try to model spin relaxation by 
the Hamiltonian Hrei = ^c|c| + R^c^^c^ , which is from 
the physical point of view dissatisfying, since it does 
not describe incoherent relaxation processes but coherent 
precession in a transverse magnetic field. This ansatz 
leads to a completely different behavior of the frequency- 
dependent current noise. Instead of a suppression of all 
spin-related effects with increasing the parameter i?, as 
expected for spin relaxation, an external field generates 
a resonance line. With increasing the field strength, this 
line just shifts to higher and higher frequencies, but does 



11 



not vanish. 



V. CONCLUSIONS 

By contacting a quantum dot to ferromagnetic leads, 
the transport characteristic through the device crucially 
depends on the quantum-dot spin. In this paper we dis- 
cussed the influence of the spin precession of the dot elec- 
tron in the tunnel-induced exchange field and an applied 
external magnetic field. While the conductance depends 
only on the time-average dot spin, the current-current 
correlation function is sensitive to its time-dependent 
evolution. 

In the zero- frequency limit, the spin precession lift the 
dynamical spin blockade, and therefore reduce the zero- 
frequency noise. At the Larmor frequency, corresponding 
to the sum of exchange and applied field, the single-spin 
precession leads to a resonance in the frequency depen- 
dent current-current correlation function. Responsible 
for the resonance is the tunnel-out process of a dot elec- 
tron to the drain lead. Due to magnetoresistance, the 
tunnel probability depend on the relative angle of dot 
spin and drain magnetization. Therefore the spin pre- 
cession leads to an oscillation of the tunnel probability, 
visible in the current-current correlation function. The 
shape of the resonance in the current-current correlation 
can either have an absorption or dispersion lineshape, 
depending on the relative angle between the lead magne- 
tizations. 

Finally, we show how to properly include spin deco- 
herence, and discuss why modelling spin relaxation by 



an external field transverse to the spin quantization axis, 
as done sometimes in the literature, is unsatisfying. 
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APPENDIX A: GENERALIZED TRANSITION 
RATES 

The generalized transition matrix W is given by the 
solution of the self energy diagrams up to linear order 
in the coupling strength F. We have choosen the quan- 
tization axis perpendicular to both lead magnetizations, 
and the x— axis symmetric with respect to the magneti- 
zations. Arranged in the matrix notation introduced in 
Sec. mTFl we get 

W|^^^=rLAL + (L^R), (Al) 

with the matrix Al given by 
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The angle (j) = 2^l = — 20r is the angle enclosed by the 
lead magnetizations. The leads are characterized by the 
Fermi functions f^{oo) and /~ = 1 — /+. For shorter 
notation we further introduced = ~ fh ~^ ^) ^ 



yi, = fi^ (e) + + [/), and the exchange field strength 
Br = \Br\/Tr, scc Eq. (|T2|) . 
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